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, Abstract. 

, This paper presents techniques developed by the LIGO Scientic Collaboration to search for 

I the stochastic gravitational-wave background using the co-located pair of LIGO interferometers 

^— ( ■ at Hanford, WA. We use correlations between interferometers and environment monitoring 

^ ' instruments, as well as time-shifts between two interferometers (described here for the first 

0^ , time) to identify correlated noise from non-gravitational sources. We veto particularly noisy 

frequency bands and assess the level of residual non-gravitational coupling that exists in the 

\l ' surviving data. 

^ ■ 

o : 

oo . 

^ I 1. Introduction 

^ ■ The LIGO (Laser Interferometer Gravitational- Wave Observatory) project has constructed 

^ ! two detectors, dubbed HI and 112, at its Hanford observatory (on the US Department of 

^ ' Energy's Hanford Site in Washington, USA) and a third, LI, at its Livingston observatory (in 



Livingston Parish, Louisiana, USA). These detectors are Fabry- Perot, power-recycled Michelson 
interferometers that detect minute changes in the differential lengths of their arms, the signatures 
of incident gravitational waves. The Fabry- Perot cavities, or arms, are 4 km long in the HI and 
LI detectors and 2 km long in H2. In November 2005, the LIGO instruments achieved their 
design sensitivities [1] and began taking data for LIGO's fifth science run, S5. The S5 observing 
effort continued through September 2007 and for some of its duration included the GEO 600 
and VIRGO instruments. This run enables searches for different types of gravitational radiation 
with unprecedented sensitivities. 

One of the LSC's (LIGO Scientific Collaboration's) major efforts is to search for an 
SGWB (stochastic gravitational- wave background), which is the gravitational analogue to the 
electromagnetic CMB (cosmic microwave background). The SGWB is predicted to contain 
contributions from unresolved populations of distant astrophysical processes involving compact 
objects as well as from early cosmological phenomena such as inflation and cosmic string 
collisions. The LSC has not yet finalized results for an SGWB search of the S5 data, but 
has published upper limits using S4 and earlier data sets ^ [Sj [4J • 



1.1. The stochastic search and H1-H2 

The SGWB is usually described in terms of the GW spectrum: 



d/ _ dpGW 
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where dpcw is the energy density of gravitational radiation contained in the frequency range / 
to / + d/ [5] , Pc is the critical energy density of the Universe, and / is frequency. 

LSC searches for the isotropic component of the SGWB deploy the cross-correlation technique 
described in Allen and Romano[5j. The technique estimates the contribution to the detector 
noise due to the SGWB by cross-correlating the detector output between a pair of detectors; 
if the detectors' intrinsic noises are uncorrelated, the common astrophysical component can be 
distinguished. We define the cross-correlation estimator: 
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Y= df Y{f) = df df 6Tif - /') slif) Uf) Qif) , (2) 
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where 5t is a finite-time approximation to the Dirac delta function, si and §2 are the Fourier 
transforms of the strain time-series of two interferometers, and Q is a filter function. In the limit 
when the detector noise is much larger than the GW signal, and assuming that the detector noise 
is stationary, Gaussian, and uncorrelated between the two interferometers, the variance of the 
estimator Y is given by: 

df aUf) « 2 X dfP,{m{f)\Q{f)\' , (3) 

where the Pi{f) are the one-sided power spectral densities (PSDs) of the two interferometers 
and T is the (single-segment) integration time. Optimization of the signal-to-noise ratio leads 
to the following form of the optimal filter [5] : 

rSi f\ AT 7(/)'S'gw(/) , n I f\ 3ifo f^Gw(/) 

Here, SQy/{f) is the strain power spectrum of the SGWB for which we are searching. Assuming 
a power-law template GW spectrum with index a, ^cwif) = ^a(//100 Hz)", the normalization 
constant M is chosen such that (Y) = Vt^T. Finally, 7(/) is the overlap reduction function. It 
captures the signal reduction due to two effects: (i) relative translation between interferometers 
introduces a frequency- and direction-dependent phase difference, and (ii) relative rotation makes 
each detector sensitive to a different polarization of incident radiation. As shown in figure [H 
the identical antenna patterns of the collocated Hanford interferometers imply 7(/) = 1, while 
for the Hl-Ll pair, the overlap reduction is significantly lower than 1 above 50 Hz. In the most 
common search case, which assumes a frequency-independent spectrum (a = 0), the theoretical 
advantage for H1-H2 over Hl-Ll is approximately a factor of ten in as shown in figured) 
This advantage could be eroded, however, by instrumental correlations between HI and H2. 

Since HI and H2 share the same vacuum enclosures and many of their optical benches 
share the same building, they are susceptible to common environmental disturbances. The 
cross-correlation induced by the environment could be greater than the noise fioor of our 
measurements, and could mimic, mask, or otherwise cause us to misestimate the amplitude 
or shape of a gravitational-wave signal. Hence, the H1-H2 pair has not been used for the SGWB 
searches so far. 

We have developed two techniques to identify non-gravitational-wave correlations between HI 
and H2 data. The IFO-PEM coherence technique was discussed in [6j. Here, IFO is shorthand 
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Figure 1. Left: Overlap reduction for the H1-H2 (black solid) and for the Hl-Ll pairs (gray 
dashed). Right: Most recent upper limits (S3 (thin solid) and S4 (thick solid) runs) on f^Q, as a 
function of a, and sensitivity predictions for S5 for the Hl-Ll (dashed) and H1-H2 (dot-dashed) 
pairs, assuming no instrumental correlations, design sensitivity, and one year of integration. 



for "interferometer" and PEM is an acronym for "physical environmental monitor". There 
are numerous sensors situated throughout the LIGO observatories providing data to the PEM 
channels; there are accelerometers, magnetometers, seismometers, voltmeters, and anemometers. 
By computing the coherence spectrum between each IFO channel and each PEM channel, one can 
find the linear environmental coupling at each frequency. The technique is sensitive to features of 
arbitrary bandwidths and resolutions, but is restricted to linear couplings and obviously cannot 
give information where PEM coverage is non-existent. 

Section [2] of this paper will describe the second, complementary handle on non-gravitational 
coupling, the time-shift technique. The time-shift technique's added value is that it is not 
hindered by incomplete knowledge of the environment; its trade-off is in its narrow bandwidth. 
Section [3] will connect the two methods to each other and to the measured quantities in the 
SGWB search and section [5] will describe how we envision applying these techniques together to 
extract the best possible SGWB measurement from the H1-H2 pair. The paper will close with 
a statement of where we stand on the H1-H2 problem. 

2. Time-shift technique 

The main idea of the time-shift technique is to perform the search described in the previous 
section after introducing an unphysical time-shift between the two interferometers' data-streams. 
That is, since the SGWB is expected to be broadband (at least hundreds of Hz wide), the 
SGWB-induced cross-correlation between the two interferometers is expected to disappear after 
a time-shift r larger than ~ 10 ms and only narrowband features, of order 1/r wide or narrower, 
are expected to be present in the time-shifted estimates of Y(f). Note that such narrowband 
features could be instrumental/environmental in nature, but they could also be genuine GW 
signals. Since narrowband GW signals are not a target of this search, our approach is to 
remove all narrow-band features from our analysis, regardless of their origin. Further note that 
a narrowband feature could also be a statistical fluctuation. In this case, it would not appear 
in the analysis with different time-shifts, and it would then have to be included in the search 
for the SGWB. While the time-shift technique is not susceptible to the "incomplete coverage" 
issues that are relevant for the IFO-PEM coherence technique, it cannot identify broadband 
instrumental correlations. Hence, the two techniques are largely complementary. 



In order to avoid correlations between estimates with different time-shifts r, the analysis 
must be performed using segment duration T < t. Since we are interested in removing ~ 1 Hz 
features and narrower, we use time-shifts of order 1 s. This also leads to the choice T = 1 s. Such 
short segment duration leads to other differences from the standard SGWB search procedure 
(where T~ 100s), described in the references [2][3][4j. In particular, the PSD estimates are 
obtained using seven neighboring Is segments on each side (instead of only one), and the data 
quality cuts are imposed on the PSDs rather than on estimates of cry, as they are less susceptible 
to downward fluctuations in the noise. We also remove the times identified by the LSC Glitch 
Group to contain instrumentally caused transients in the data. The method was tested using the 
S5 data between November 2005 and April 2006, using the calibrated strain data, a 40-500 Hz 
frequency range, a sixth-order butterworth high-pass filter at 32 Hz to remove low-frequency 
detector noise, and 50% overlapping Hann windows. In the study, we never estimated Y{f) 
without an unphysical time-shift, to avoid potentially biasing the resulting data-quality cuts. 

We define the frequency-dependent signal-to-noise ratio, as a function of time-shift r, properly 
accounting for the frequency bin width df: 

o-(/;r) Vd/ 

Using Is segments, we performed ibl,ib2s time-shifts. We then attempted three ways 
of combining these four measurements to estimate the the SNR at zero-shift, SNRts(/;0): 
Gaussian fit to the four measurements, maximum over the four measurements, and average of 
ibis measurements. These three approaches led to very similar results, so in the following we 
conservatively use the maximum over the three. 

3. Comparison of the two techniques 

The PEM-IFO coherence method estimates the environmental contribution to the coherence 
between the two interferometers 7pem(/)- Coherence 7(/) is defined by: 

-Yn- /i2(/) _ Pi2(.f)Q{f)df 1 _ SNR(/) 

VPiimif) VA(/)P2(/)Q2(/)d//T vTdf VTdf ' 

where Puif) = (si(/)s2(/)) is the cross-power spectral density and T is the total observation 
time, so Tdf is the total number of segments. Hence, we can estimate the environmental 
contribution to the SNR of the stochastic search: 



SNRpEM (/) = 7PEM (/) \/Td7. (7) 

This allows us to directly compare the two techniques, keeping in mind that SNRts is defined 
using |y(/)|. Moreover, we can define the environmental contribution to the point estimate, 
where again, the PEM subscript indicates that the quantity is estimated using the PEM-IFO 
coherence method: 

n yPEM(/) 7PEm(/)^(/) 
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4. Proposed search algorithm 
4.1. Veto 

The first obvious application of the described techniques is to veto frequencies in which HI and 
H2 exhibit strong non-gravitational coupling. Figure [2] compares the SNR estimates using the 
two techniques. Note that the two techniques generally agree well in identifying contaminated 
frequency bands. These bands will be vetoed in the SGWB search. It is not completely self- 
evident precisely how this veto should be defined, since the IFO channels are coupled to some 
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Figure 2. Comparison of the PEM-IFO coherence and time-shift techniques for estimating the 
contaminated frequency bands. The 60 Hz harmonics have ah'eady been vetoed. 



extent at all frequencies. One approach is to veto frequency bins in which both methods' SNRs 
exceed 2. This choice is somewhat relaxed, and it relies on the assumption that we can estimate 
r^PEM in the remaining frequency bands that pass the veto, as discussed below. 



Estimate the residual environmental contribution 
Using equation [U we can estimate the environmental contribution to the estimate of the SGWB 
amplitude r^pEMi in the frequency bands that pass the veto described in section 14.11 This 
estimate could be potentially subtracted from the overall SGWB amplitude estimate, thereby 
producing a better estimate of the gravitational-wave contribution. 

In this case, it is important to understand the uncertainty in the estimate of ilpEM- The 
statistical uncertainty in flpEM comes from the statistical uncertainties on 7pem and on ay, 
both of which are expected to be negligible compared to the theoretical uncertainty on the 
overall amplitude estimate an = cry/T. 



4-3. Estimate systematics 

The systematic error on $7pem could potentially be large. A few major sources of uncertainty 
are apparent: the approximation made in the definition of 7pem from the individual IFO-PEM 
coherences [6j and differences in the details of the data manipulation between the IFO-PEM 
coherence and SGWB search calculations. One approach for estimating this systematic error 
is to compare the f^pEM estimate with the SGWB estimate in a frequency band that is known 
to be dominated by environmental correlations. As the environment dominates, fipEM should 
explain the entire SGWB estimate of in this band. The mismatch would give us an estimate 
on the systematic error of our IFO-PEM coherence measurement. 

We have performed such a calculation. Our preliminary results indicate that the IFO-PEM 
and the SGWB estimates agree within 50% in the contaminated frequency bands. However, we 
believe we can significantly improve this agreement; the differences in the data manipulation in 
the two techniques (segment duration, frequency resolution, PSD weighting, etc) were significant 
and we expect to minimize them in future iterations. 

Another indication of the systematic error could be obtained from hardware injections. 
At semi-regular intervals in the LIGO science runs, we injected simulated SGWB signals by 
physically shaking the optics; the SGWB search on these segments tests how accurately we can 
recover the injected signal's amplitude. If we instead treat the injection data stream as a PEM 



channel, we can verify that we can subtract out the injection and instead produce an estimate 
consistent with noise. A signal detection would imply an imperfect subtraction. 

4.4. SGWB search 

Having identified the most egregiously correlated frequency bands, and estimated the 

environmental contribution fipEMj we can proceed with the usual SGWB search algorithm. 
This includes a number of diagnostics on the quality of data, such as tests of Gaussianity, long- 
term stationarity, etc. Potential problems would become apparent here, such as omissions in the 
frequency and time vetoes. Any problems would be investigated and used as input in another 
iteration of the veto and environment contribution estimation steps. The very last step is to 
produce the SGWB estimate, and potentially subtract the ripEM estimate. 

5. Conclusions 

We have presented the time-slide technique, by which all narrowband, correlated noise between 
two interferometers can be identified. It provides complementary coverage of the space of possible 
non-gravitational couplings between the interferometers, alongside the IFO-PEM coherence 
technique. It is important to note that these techniques leave some blind spots: they leave 
out the possibility of broad-band environmental contributions that are unmonitored in the 
observatory's environment, and they misestimate their contributions if the couplings are non- 
linear. This includes effects such as seismic upconversion, a non-linear process by which low- 
frequency seismic activity excites higher- frequency vibrations in the instrument, and stray light 
reflections in the beam tubes that introduce cross-talk between the interferometers. 
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